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Hierarchic  elements  which  implement  the  p-version  efficiently  are  used  to¬ 
gether  with  precomputed  arrays  of  elemental  stiffness  matrices. 

Some  important  research  results  are: 

1.  In  polygonal  regions,  under  conditions  usually  satisfied  in  practice,  if 
the  h-version  converges  with  order  of  error  in  energy  norm  OCl/N01)  then 
the  p-version  converges  with  order  of  error  in  energy  norm  0(1/N^a) , 
where  N  is  the  number  of  degrees  of  freedom.  This  applies  to  planar 
problems  which  require  either  CO  or  Cl  global  continuity. 

2.  Hierarchic  CO  elements  for  triangles  or  rectangles  have  been  developed 

and  implemented  for  p^2 .  Hierarchic  Cl  triangular  elements  have  been 
developed  and  implemented  for  p_>5.  Hierarchic  CO  solid  elements  have 
been  developed  for  p>2  for  bricks,  tetrahedra,  triangular  prisms,  and 
rectangular  pyramids .  "  . 

3.  Modified  Bernstein  polynomials  have  been  constructed  over  triangles  which 

provide  a  smooth  approximation  to  functions  in  H^.  These  polynomials  are 
used  to  prove  that  the  p-version  of  the  finite  element  method  converges 
in  the  Cl  case  (plate  bending) .  ■ 

4.  Explicit  mappings  have  been  constructed  which  map  triangles  with  one 
curved  side  into  the  standard  triangle.  These  mappings  have  the  property 
that  the  resulting  elemental  stiffness  matrices  can  be  integrated  in 
closed  form.  The  need  for  numerical  quadrature,  which  may  be  inefficient 
in  the  p-version,  is  thereby  obviated.  The  curved  side  can  be  either 
parabolic  or  elliptic.  The  newly  constructed  mappings  permit  the  p-ver¬ 
sion  to  be  applied  to  domains  with  curved  boundaries  of  specified  shapes. 

5.  The  p-version  analysis  of  stresses  in  a  rhombic  plate  was  completed. 

The  results  compare  very  favorably  with  those  of  other  codes. 
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1.  INTRODUCTION 


Two  approaches  to  finite  element  analysis  are  now  widely  recognized  in  the 
engineering  and  mathematical  communities.  In  both  approaches  the  domain  0  Is 
divided  Into  simple  convex  subdomains  (usually  triangles  or  rectangles  in  two 
dimensions,  and  tetrahedra  or  bricks  in  three  dimensions)  and  over  each  sub- 
domain  the  unknown  (displacement  field)  is  approximated  by  a  (local)  basis 
function  (usually  a  polynomial  of  degree  p) .  Basis  functions  are  required  to 
join  continuously  at  boundaries  of  the  subdomains  in  the  case  of  planar  or  3 
dimensional  elasticity,  or  smoothly  in  the  case  of  plate  bending.  The  differ¬ 
ence  between  the  two  approaches  lies  In  the  manner  In  which  convergence  is 
achieved.  These  two  approaches  are: 

1.  The  h-version  of  the  finite  element  method.  In  this  approach  the  degree  p 
of  the  approximating  polynomial  Is  kept  fixed,  usually  at  some  low  number 
such  as  2  or  3.  Convergence  is  achieved  by  allowing  h,  the  maximum  diameter 
of  the  convex  subdomains,  to  go  to  zero.  Estimates  for  the  error  in  energy 
have  long  been  kr  ;. . [1,  2,  3].  In  all  of  these  estimates  p  is  assumed  to 
be  fixed  and  the  error  estimate  is  asymptotic  in  h,  as  h  goes  to  zero. 

2.  The  p-version  of  the  finite  element  method.  In  this  approach  the  subdivi¬ 
sion  of  the  domain  (2  is  kept  fixed  but  p  is  allowed  to  increase  until  a 
desired  accuracy  is  attained.  The  p-version  is  reminiscent  of  the  Rltz 
method  for  solving  partial  differential  equations  but  with  a  crucial  distinc 
tlon  between  the  two  methods.  In  the  Rltz  method  a  single  polynomial 
approximation  is  used  over  the  entire  domain  ft  (ft.  In  general.  Is  not  con¬ 
vex)  .  In  the  p-version  of  the  finite  element  method  polynomials  are  used 

as  approximations  over  convex  subdomains.  This  critical  difference  gives 
the  p-version  a  much  more  rapid  rate  of  convergence  than  either  the  Ritz 
method  or  the  h-version  as  will  be  explained  later. 


The  p-version  of  the  fiuite  element  method  requires  families  of  polyno¬ 
mials  of  arbitrary  degree  p  defined  over  different  geometric  shapes.  Polyno¬ 
mials  defined  over  neighboring  elements  join  either  continuously  (are  in  C^) 
for  planar  or  three  dimensional  elasticity,  and  smoothly  (are  in  C^)  for  plate 
bending.  In  order  to  implement  the  p-verslon  efficiently  on  the  computer, 
these  families  should  have  the  property  that  computations  performed  for  an 
approximation  of  degree  p  are  re-usable  for  computations  performed  for  the 
next  approximation  of  degree  p  +  1.  We  call  families  possessing  this  property 
hierarchic  families  of  finite  elements. 

The  h-version  of  the  finite  element  method  has  been  the  subject  of  inten¬ 
sive  study  since  the  early  1950's  and  perhaps  even  earlier.  Study  of  the 
p-version  of  the  finite  element  method,  on  the  other  hand,  began  at  Washington 
University  in  St.  Louis  in  the  early  1970's.  Research  in  the  p-version  (formerly 
called  The  Constraint  Method)  has  been  supported  in  part  of  the  Air  Force  Office 
of  Scientific  Research  since  1976.  In  this  report  we  review  some  of  the  more 
important,  accomplishments  during  this  funding  period,  and  describe  the  current 
status  of  the  p-verslon  of  the  finite  element  method. 

2.  RESEARCH  ACCOMPLISHMENTS 
2.1.  Rate  of  convergence  of  the  p-version 

Extensive  computational  experiments  have  provided  empirical  evidence 
that  the  rate  of  convergence  in  the  p-version  is  significantly  higher  than 
in  the  h-verslon.  For  example,  a  shell  problem  consisting  of  a  circular 
cylindrical  shell  with  symmetrically  located  cutouts  was  subjected  to  a 


uniform  axial  end  shortening  of  known  amount  as  shown  below. 


I 

(ITU* 


koexMtto  rcsr  moun  N*t 


rUHO  V  HMMCTKr 


xcnoN  •■ 


The  boundary  conditions  at  the  ends  of  the  shell  are 
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u  ■  constant  ■  0.2  x  10  Inches 


The  solution  to  this  problem  using  the  p-version  of  the  finite  element  was 
compared  to  the  solution  by  other  computer  programs  in  [4].  The  table  shown 
below  gives  the  number  of  degrees  of  freedom  needed  by  the  different  computer 
codes  to  attain  the  same  degree  of  accuracy. 

It  is  clear  that  the  p-version  requires  far  fewer  degrees  of  freedom  than  the 
next  best  code,  TRISHL  which  was  specifically  altered  to  solve  this  problem. 


COMPUTES 

CODE 


SOURCE 


NUMBER  OF  DEGREES  OF 
ELEMENTS.  FREEDOM 


SELL  9 

SUIT  SIN  UAL  ATOMIC, 

INC. 

475 

2457  . 

STAGS 

LOCKHEED  MISSILES  AND 
SPACE  CORN. 

342 

1437 

REXBAT 

LOCKHEED  MISSILES  AMO 
SPACE  COAT. 

241 

APPROX  1125 

TRISHL 

NATIONAL  AERONAUTICAL 
ESTABLISHMENT,  CANADA. 

100 

637 

CDflET-X 

WASHINGTON  UNIVERSITY' 

10 

500 

(6-6-7-caseJ 

COMPARISON  OP  SOM  COMPUTATIONAL  PARAMETERS 


In  other  sample  problems  the  reduction  of  the  number  of  degrees  of  freedom 
used  by  the  p-version  Is  even  more  striking.  In  [5],  the  p-version  is  applied 
to  problems  in  elastic  fracture  mechanics  with  excellent  results  and  in  [6] 
the  p- version  is  used  to  analyze  an  edge-cracked  panel  and  a  parabolically 
loaded  panel.  In  Figures  1  and  2  the  respective  triangulations  are  shown  and 
in  Figures  3  and  4,  the  error  in  strain  energy  is  plotted  against  the  recipro¬ 
cal  of  the  number  of  degrees  of  freedom. 

The  following  two  theorems,  which  were  recently  proved,  provide  a 
rigorous  explanation  for  the  efficiency  of  the  p-version.  In  both  theorems 
ft  Is  a  bounded  polygonal  domain  in  the  plane.  In  Theorem  1,  a  model  problem 
for  the  case  is  considered,  and  In  Theorem  2,  a  model  problem  for  the 
case  is  considered.  In  both  cases,  the  problems  are  singularity  problems, 
that  is  the  smoothness  of  the  solutions  are  governed  by  the  local  behavior  at 
the  vertices  A^  of  the  polygons.  Suppose  that  is  the  angle  at  vertex 
and  that  polar  coordinates  (r^,  are  used  at  A^. 
neighborhood  of  A^  is  of  the  form 


The  solution  In  the 
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pi(ri>  W 

Yi 

Pi(ri)  -  r±  g±(|log  rj),  0^^)  is  very  smooth, 
depends  upon  a^. 

The  domain  £2  is  assumed  to  be  triangulated  in  such  a  way  that  coincides 
with  a  vertex  of  a  triangle.  The  exact  assumptions  on  g^  and  6^  are  given 
In  [6] . 

Theorem  1  (model  problem  for  Case) 


-Au  +  u  “  f  in  n 


u  ■  0  on  30 


Let  u  be  the  solution  to  (1)  in  the  weak  sense,  and  let  u  be  the  finite 

P 

element  approximation  to  u,  using  polynomials  of  degree  p  with  the  trlangu- 
latlon  S  fixed  (l.e.  up  is  the  solution  to  (1)  using  the  p-version  of  the 
finite  element  method).  If  u  e  H  (Q),  with  k  >  1,  then 

I  k  -  up I  I2>q  1  c  I M  lk»  w  ■  0*-i»  2^) ,  Yi  “  ^  (2) 

where  e  >  0  is  arbitrary. 

Theorem  2  (model  problem  for  case) 
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w  «  - —  *  0  on  an  (clamped  edge) 
dn 

Let  w  be  the  solution  to  (2)  in  the  weak  sense,  and  let  w  be  the  solution  to 

P 

(3)  using  the  p-version  of  the  finite  element  method.  If  w  e  H  (n)  with 
k  >  2,  then 


6- 


|  |v  -  wpM2ffl  1  c  p-t,+€  I  lwl  lk»  W  “  min(k-2,  2(y±  -  D)  (4) 

Y^(a^)  ■  smallest  positive  root  of  the  equation 
sln2^  -  l)ai  -  (y±  -  l)2  sin2  a±  -  0 


The  error  estimates  In  (2)  and  (4)  can  be  compared  with  analogous  estl 
mates  for  u^  and  w^,  the  solutions  by  the  h-version  of  the  finite  element 
method.  Assume,  for  convenience  that  in  (2)  k  -  1  >_  2y^»  and  in  (4) 
k  -  2  >_  2  -  1),  that  is  convergence  is  determined  by  the  nature  of  the 

singularities  at  corners.  Then  the  analogous  estimates  are 


i  1“  -  “hi  lifn  -  ch  YHullk»  Y  -  min  yx  (2’) 

•  lW  *  Wh>  l2,n  -  Ch<Y“1>  I  lwi  Ik*  Y  -  min  y±  •  (4’) 

If  N  is  the  number  of  degrees  of  freedom  then  p  -  N1^2,  h  -  N-1/2  and 

1 1"  -  “pHi.a  '  0<k'T)>  1 1"  '  "hi  li,a  ’  0(n'y/2) 

1 1»  -  »pl l2>1)  -  o(»-(T-l)),  1 1*  -  *hl l2>fl  -  oar1™*-”) 

Therefore,  if  the  criterion  used  to  compare  methods  is  the  number  of  degrees 
of  freedom  required  to  achieve  a  given  error  in  energy,  then  the  rate  of  con¬ 
vergence  of  the  p-version  is  twice  that  of  the  h-version.  This  result  provides 
a  rigorous  proof  for  the  extensive  computational  evidence  that  has  been  gathered, 
and  explains  in  part  the  efficiency  of  the  p-version. 

2.2.  Hierarchic  Families  of  Solid  Finite  Elements 

In  order  to  Implement  the  p-version  efficiently,  families  of  finite  elements 
are  needed  with  the  hierarchic  property:  computations  performed  for  an  approxi¬ 
mation  of  order  p  should  be  re-usable  when  raising  the  order  to  p  +  1.  More 


specif icially,  Che  stiffness  matrix  corresponding  to  the  polynomial  approxi¬ 
mation  of  degree  p  should  be  a  submatrix  of  the  polynomial  approximation  of 
degree  p  +  1.  In  terms  of  basis  functions,  this  implies  that  the  basis  func¬ 
tions  for  a  pth  order  approximation  should  be  a  subset  of  the  basis  functions 
for  a  (p  +  l)st  order  approximation. 

Hierarchic  families  for  triangles  both  in  the  case  and  in  the  case 
are  described  in  detail  in  [7,  8,  9,  10,  11].  We  now  describe  hierarchic 
families  of  polynomials  for  various  three  dimensional  shapes.  All  of  these 
families  are  globally  in  C^. 

tetrahedron.  A  hierarchic  family  for  the  tetrahedron  can  be  constructed  from 
the  hierarchic  family  for  the  triangle  by  using  natural  coordinates.  The 
linear  element  is: 


vertex 

nodes 


nodal 

variable 

shape  function 

u(l) 

L1 

u(2) 

L2 

u(3) 

L3 

u(4) 

L4 

4 


The  quadratic  element  has  the  same  nodal  variables  and  shape  functions  as  the 


linear  element  and  the  additional  nodes: 


nodal 

variables 

shape  functions 

uss  <U> 

“  2  L1  L2 

<“> 

“  2  L2  L3 

u88  <31> 

-lllL3 

uss  <14> 

- 1  li  H 

uss  <34> 

-  1  L3  L4 

uss  <24> 

*  7  L2  l4 

The  hierarchic  family  for  a  tetrahedron  up  to  the  element  of  degree  4 
la  shown  In  Figure  5.  Boundary  nodal  variables  correspond  to  vertices,  edges, 
and  faces.  Additional  Internal  nodal  variables  are  introducted  to  provide  a 
basis  for  a  complete  polynomial  of  degree  p.  Complete  details  are  given  in  [12] 
brick  A  hierarchic  family  for  the  brick  can  be  constructed  from  the 
hierarchic  family  for  the  rectangle.  For  convenience,  we  consider  only  the 
square  and  from  it  we  construct  a  hierarchic  family  for  the  cube.  The  hierar¬ 
chic  family  for  the  square  Is  given  as  follows. 
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HIERARCHIC  QUAORATtC 
C°  SQUARE  CUWT 


satisfies 


Q^+l)  -  0  Q*15  (0)  -  0  i  -  2,  ...,  j  -  1 

<0)  -  1. 

Now,  consider  the  square  of  side  2  as  shown  in  the  figure.  Basis  functions 
corresponding  to  the  nodal  variables  u(i)  i  ■  1,  2,  3,4  are 


"uCl)  “  J  <l  '  «>«  -  n)  \(3)  ■  7  a  +  oa  +  n) 

MuC2)  ■  T  a  +  5)(1  -  ">  Hu(»)  -ft1-  »<l  +  ") 


(5) 


and  it  is  easily  seen  that  these  basis  functions  span  the  same  space  as 
l,€,n,€n  i.e.  they  contain  the  complete  linear  polynomial.  Also  these  nodal 
variables  enforce  C°  continuity  across  sides.  Now  denoting  by  (ij)  the  mid¬ 
point  of  side  ij,  basis  functions  corresponding  to  the  nodal  variables  11^(12), 

unn<23)’  uU(34)’  unn(41)  are 

N  (12)  -  Q?(5)  (1  -  n)  N  (34)  -  Q2(-5)(l  +  n) 

(6) 


Nu  (23)  -  (1  +  C)Q2(n)  Nu  (41)  -  (1  -  5)Q2(-n) 

nn  nn 

and  these  basis  functions  added  to  the  ones  in  (5)  span  the  same  space  as 
2  2  2  2 

l»5»n,5  ,5n,n  ,5  n,5n  ,  i.e.  they  contain  a  complete  quadratic  and  they  enforce 
continuity  along  sides.  The  basis  functions  (5)  are  taken  for  the  hierar¬ 
chic  rectangular  C3  linear  element,  and  those  in  (5)  and  (6)  for  the  quadratic 
element . 

For  j  _>  3,  we  have 


Nu  j  (12)  -  QjUMI  -  n) 


N  . (34) 


Qj(-C)(l  +  n) 


N  . (23)  -  (1  +  5)Q,(n) 


N  ,(41)  -  (1  -  OQ.(-n) 


as  Che  basis  function  for  jth  order  tangential  derivatives  at  midsides,  and 
for  j  _>  4  we  add  the  internal  modes 

(1  -  C2)  (1  -  n2)  J-4~V  i  -  0 . j  -  4. 


These  basis  functions  span  the  same  space  as  a  complete  polynomial  of  degree 
and  the  two  monomials  of  degree  j+1,  5^n,  5n^ . 

Thus,  the  hierarchic  square  element  of  degree  p  >_  2  has  (1/2) (p+1) (p+2)+2 
basis  functions,  two  more  than  the  dimension  of  the  complete  polynomial  of 
degree  p.  The  two  extra  uerms  correspond  to  5Pn,  5nP.  By  scaling  the  sides 
of  the  square  the  elements  are  easily  transformed  into  rectangular  ones. 

The  hierarchic  family  for  the  brick  can  now  be  Constructed .  Linear  element: 
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The  quadratic  element  adds  the  shape  functions  as  follows: 
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Higher  degree  elements  add  edge  modes,  face  modes  and  internal  modes.  More 
details  are  given  In  [12,  13]. 

triangular  prism  A  hierarchic  family  for  the  triangular  prism  can  be 
constructed  by  combining  the  hierarchic  families  for  the  triangle  (in  natural 
coordinates)  and  for  the  square  (In  rectangular  coordinates) 
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The  quadratic  element  adds  the  following  shape  functions 
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Higher  degree  elements  add  edge  modes,  face  modes  and  internal  modes.  Further 
details  are  given  in  [12,  13]. 

square  pyramid.  It  can  be  shown  [13]  that  no  hierarchic  family  exists 
which  consists  of  polynomials  alone.  However  it  is  possible  to  supplement 
certain  rational  functions  in  such  a  way  that  the  element  of  degree  p  contains 
a  complete  polynomial  of  degree  p  and  additional  rational 


functions.  Furthermore,  because  of  the  special  form  of 
these  rational  functions,  integration  of  all  shape  func¬ 
tions  which  appear  in  the  elemental  stiffness  matrix  can 
be  performed  in  closed  form.  No  numerical  quadrature 
is  required.  The  linear  element  is: 
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Each  of  the  vertex  modes  is  a  linear  combination  of  1,  5,  n,  ?,  . 

Higher  degree  elements  again  add  edge  modes ,  face  modes  and  Internal  modes. 
Further  details  are  given  in  [13]. 

Combinations 

Hierarchic  families  of  the  different  shapes  shown  here  have  been  con¬ 
structed  so  that  they  join  together  continuously.  Thus  geometries  such  as  the 
ones  shown  below  can  be  modeled  using  as  few  elements  as  possible.  In  this 
way  the  subdivision  of  the  polyhedral  domain  ft  can  be  made  very  coarse  and 
accuracy  can  be  obtained  by  Increasing  p,  without  using  added  degrees  of 
freedom  to  describe  the  geometry.  As  has  been  shown  (at  least  for  two  dimen¬ 
sional  problems)  the  rate  of  convergence  in  the  p-version  is  twice  that 
obtained  for  the  h-version.  Therefore  subdividing 


in  this  way  leads  to  computational  efficiency,  since  it  makes  maximum  use  of 
the  p-version.  it  also  requires  significantly  simpler  input. 

2.3  P-Version  Analysis  of  a  Rhombic  Plate 


The  problem  of  the  stress  analysis  of  a  rhombic  plate  is  a  difficult  one 
which  has  attracted  much  attention  [14,  15,  16,  17].  We  consider  the  rhombic 
plate  shown  in  Figure  6,  which  is  simply  supported,  uniformly  loaded  and  whose 
acute  angles  are  30  degrees.  A  very  strong  singularity  occurs  at  the  obtuse 
vertex  when  plate  bending  theory  is  applied  to  this  problem,  and  this  singularity 
leads  to  y  -  0.4  in  Theorem  2,  Equation  (4).  As  a  result  many  finite  element 


4 


models  of  Che  rhombic  plate  either  fall  to  converge  or  converge  very  slowly. 
The  hierarchic  family  used  in  the  p-version  contains  corrective  rational 
functions.  These  rational  functions  absorb  the  singularity  very  well  and  lead 
to  the  best  results  thus  far  obtained  for  this  problem.  Shown  below  are  the 
values  for  the  deflection  at  the  central  point,  we,  using  the  p-version 
analysis,  with  p  ■  5  and  with  different  numbers  of  elements.  These  results 
are  compared  with  results  obtained  by  other  methods 
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other  codes 


no.  of  degrees  of  freedom 

%  error 

Sander  [14] 

>  1000 

24 

Argyris  [ 15 ] 

^  1300 

2 

Basu,  Szabo  et  al  [16] 

~  130 

2 

Basu,  Szabo  et  al  [16] 

:  180 

i 

Morley  [17]  exact  value  (analytic  expression)  =  4.08 


COMET-X  is  an  experimental  computer  code  which  implements  the  p -vers ion 
of  the  finite  element  method  by  using  the  hierarchic  families  which  have  been 
constructed.  COMET-X  is  maintained  by  the  Center  for  Computational  Mechanics 
at  Washington  University  and  it  is  the  only  computer  code  In  existence  which 
Implements  the  p-version.  COMET-X  can  be  used  as  a  code  to  Implement  the  h- 
verslon  as  well  simply  by  fixing  the  polynomial  order  p  and  refining  the  mesh. 

COMET-X  currently  has  the  following  capabilities: 

A.  Element  types :  Stiffeners,  triangular  elements,  triangular  elements 
with  one  side  curved,  rectangular  elements,  solid  elements  of  the  shapes 
described  earlier. 

B.  Types  of  Analysis:  Laplace  and  Poisson  equations,  plane  elasticity, 
elastic  plate  bending,  stiffened  elastic  plates,  3 -dimensional  elasticity, 
temperature  distribution  in  3 -dimensions. 

C.  Special  Capabilities:  non-uniform  p-distrlbution,  elastic  fracture 
mechanics  computations  in  two  dimensions,  nearly  incompressible  solids, 
linear  boundary  layer  problems. 

D.  Pre-  and  Post  processing  capabilities  Including  graphics  and  visual  displays. 
The  capabilities  and  usage  of  COMET-X  are  described  in  detail  in  [ 1 6 ] . 


As  an  example  of  a  problem  analyzed  by  the  p-version  using  hierarchic 
solid  finite  elements,  [18],  consider  the  finite  element  model  of  a  gear 
casing,  shown  below.  The  model  on  the  right  was  used  in  COMET-X,  the  one  on 
the  left  was  used  by  NAS TRAN.  Stress  contours  shown  in  Figure  7  were  obtained 
from  COMET-X. 


COMET-X  NASTRAN 


Finite  Element  Meshes  for  a  Gear  Casing 
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Figure  6.  The  Rhombic  Plate  Problem 
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abstract 


With  one  exception,  all  finite  element  software  systems  have  element 
libraries  in  which  the  approximation  properties  of  elements  are  frozen. 

The  user  controls  only  the  number  and  distribution  of  finite  elements. 

The  exception  is  an  experimental  software  system,  developed  at  Washington 
University.  This  system,  called  COMET-X,  employs  conforming  elements 
based  on  complete  polynomials  of  arbitrary  order.  The  elements  are  hierar¬ 
chic,  i.e.  the  stiffness  matrix  of  each  element  is  embedded  in  the  stiff¬ 
ness  matrices  of  all  higher  order  elements  of  the  same  kind.  The  user  con¬ 
trols  not  only  the  number  and  distribution  of  finite  elements  but  their 
approximation  properties  as  well.  Thus  convergence  can  be  achieved  on 
fixed  mesh.  This  provides  for  very  efficient  and  highly  accurate  approxi¬ 
mation  and  a  new  method  for  computing  stress  intensity  factors  in  linear 
elastic  fracture  mechanics.  The  theoretical  developments  are  outlined, 
numerical  examples  are  given  and  the  concept  of  an  advanced  self-adaptive 
finite  element  software  system  is  presented. 


2.  "The  Constraint  Method  for  Finite  Element  Stress  Analysis",  by  I.  N. 

Katz,  presented  at  the  National  Bureau  of  Standards,  Applied  Mathe¬ 
matics  Dlvlson  on  October  19,  1977. 

abstract 

In  conventional  approaches  to  finite  element  stress  analysis  accuracy 
is  obtained  by  fixing  the  degree  p  of  the  approximating  polynomial  and  by 
allowing  the  mn-jHmim  diameter  h  of  elements  in  the  triangulation  to  approach 
zero.  An  alternate  approach  is  to  fix  the  trlangulatlon  and  to  increase  the 
degrees  of  approximating  polynomials  in  those  elements  where  more  accuracy 
is  required.  In  order  to  Implement  the  second  approach  efficiently  it  is 
necessary  to  have  a  family  of  finite  elements  of  arbitrary  polynomial  degree 
p  with  the  property  that  as  much  information  as  possible  can  be  retained  from 
the  pth  degree  approximation  when  computing  the  (p+1) st -degree  approximation. 
Such  a  HIERARCHIC  family  has  been  formulated  with  p  2  for  problems  in  plane 
stress  analysis  and  with  p  >_  5  for  problems  in  plate  bending.  The  family  is 
described  and  numerical  examples  are  presented  which  illustrate  the  efficiency 
of  the  new  method. 


3.  "The  p-version  of  the  Finite  Element  Method",  by  I.  N.  Katz  and  B.  A. 

Szabo,  presented  at  Air-Force  Flight  Dynamics  Laboratory  Wright- 

Patterson  Air  Force  Baae  on  April  23,  1981. 

abstract 

The  theoretical  basis  of  the  p-version  of  the  finite  element  method  has 
been  established  only  quite  recently.  Nevertheless,  the  p-version  is  already 
seen  to  be  the  most  promising  approach  for  Implementing  adaptivity  in  practi¬ 
cal  computations.  The  main  theorems  establishing  asymptotic  rates  of  conver¬ 
gence  for  the  p-version,  some  aspects  of  the  algorithmic  structure  of  p-verslon 
computer  codes,  numerical  experience  and  a  posteriori  error  estimation  will  be 
discussed  from  the  mathematical  and  engineering  points  of  view. 


